ADA053526 


0 


SSS-R-77-3057 


THE  DEPENDENCE  OF  BODY  WAVE  MAGNITUDE  ON  YIELD 
FOR  UNDERGROUND  EXPLOSIONS  IN  SALT 

Topical  Report 

T.  C.  Bache 
J . T . Cherry 
B.  F.  Mason 


Sponsored  by 

Advanced  Research  Projects  Agency 
ARPA  Order  No.  2551 

This  research  was  supported  by  the  Advanced 
Research  Projects  Agency  of  the  Department  of 
Defense  and  was  monitored  by  AFTAC/VSC,  Patrick 
AFB,  FL  32925,  under  Contract  No.  F08606-76-C- 
0041. 

The  views  and  conclusions  contained  in  this  docu- 
ment are  those  of  the  authors  and  should  not  be 
interpreted  as  necessarily  representing  the  of- 
ficial policies,  either  expressed  or  implied,  of 
the  Advanced  Research  Projects  Agency,  the  Air 
Force  Technical  Applications  Center,  or  the  U.S. 
Government . 

Approved  for  Public  Release,  Distribution  Unlimited. 


November  1976 


P.  0.  BOX  1620,  LA  JOLLA,  CALIFORNIA  92038,  TELEPHONE  1714)  463-0060 


AFTAC  Project  Authorization  No.  VEIjA/T/7712/B/ETR 
Program  Code  No.  6H189 

Effective  Date  of  Contract:  October  1,  1976 

Contract  Expiration  Date:  September  30,  1977 

Amount  of  Contract:  $374,397 

Contract  No.  F08606-76-C-0041 

Principal  Investigator  and  Phone  No. 

Dr.  Thomas  C.  Bache,  (714)  453-0060,  Ext.  337 

Project  Scientist  and  Phone  No. 

Dr.  Ralph  W.  Alewine,  III,  (202)  325-8484 


UNCLASSIFIED 

SCCuKITV  CLASSIFICATION  OF  THIS  PACE  Oat*  Enfrtd) 


REPORT  DOCUMENTATION  PAGE 


5S5--  ii-77-3<^5-7 


2.  GOVT  ACCESSION  N 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RCCIPICNT'S  catalog  NUMSER 


I li:|ra*l.(rfT 


NIZATION  name  and  AODPESS 

Systems,  Science  emd  Software  ^ 

P.  0.  Box  1620 

La  Jolla,  California  92038 


IV  CONTNOLLING  OFFICE  NAME  AND  AOOPESS 

VELA  Seismological  Center 
312  Montgomery  Street 
Alexandria,  Virginia  22314 


I 


MONtTORlNG  AGENCY  NAME  S AODRESS^I/  dUfmttxt  Inm  Controltlng  OHic 


I 


rhis  reporf) 


Unclassified 


ISa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


1«.  distribution  STATEMENT  (ol  thit  Raport; 


Approved  for  Public  Release,  Distribution  Unlimited. 


17.  OfSTRlEuTlON  STATEMENT  (of  III#  mb9tr»et  Mfrmd  in  Block  20,  if  difformnt  from  Roport) 


It-  KEY  WORDS  (Continuo  on  rovorco  Mido  if  nocottorr  ^d  Idontity  by  block  number) 

Nuclear  Explosion  Seismology 

Body  Wave  Magnitude  Predictions  lO 

xplosion  Source  Theory 


20.  ABSTN^T  fConflmi*  on  rovoroo  aldo  If  nocosaoiy  ontf  Idontity  by  block  number)  J 

The  dependence  of  body  wave  magnitude  (mb)  on  yield  for 
explosions  in  salt  is  addressed.  The  objective  is  to  compare 
Eurasian  explosions  in  salt  to  hypothetical  explosions  in  salt 
at  Nevada  Test  Site  (NTS)  and  to  actual  NTS  explosions  in 
granite.  The  coupling  of  the  explosion  energy  into  elastic 
waves  is  computed  with  a finite  difference  program.  Four  cal- 
culations were  done  in  which  all  material  properties  were  fixed-.^ 


DD  1473  edition  of  I NOV  (S  IS  oesolete 


3$^  7 


UNCLASSIFIED 

security  classification  of  This  page  flWi»«  Dmtm  Enlttad) 


R-3057 


TABLE  OF  CONTENTS 


I.  INTRODUCTION  

II.  THE  EQUIVALENT  ELASTIC  SOURCE  FOR  EXPLOSIONS 

IN  SALT  

III.  THE  DEPENDENCE  OF  m.  ON  YIELD  FOR  EXPLOSIONS 

IN  Si\LT 

3.1  COMPUTATIONAL  METHOD  

3.2  EARTH  MODELS  

3.3  SYNTHETIC  SEISMOGRAMS  

3.4  THE  DEPENDENCE  OF  nij^  ON  YIELD 

3.5  COMPARISON  TO  OBSERVATIONS  

3.6  COMPARISON  TO  HYPOTHETICAL  SALT  EVENTS 

AT  NTS  

3.7  COMPARISON  TO  EVENTS  IN  GRANITE  . . . 1' 

3.8  SALMON i . . 

REFERENCES  


Page 

1 

3 

8 

8 

9 

11 

13 

17 

17 
21 
27 
35  , 


R-3057 


I.  INTRODUCTION 

In  this  report  we  address  the  dependence  of  body  wave 
magnitude  (m^)  on  yield  for  explosions  in  salt.  Our  objective 
is  to  compare  salt  events  in  Eurasia  to  similar,  though  hypo- 
thetical, events  at  NTS  and  to  granite  events  at  NTS,  Our 
study  uses  deterministic  computational  methods  and  our  re- 
sults are  based  on  time  domain  measurements  from  synthetic 
seismograms.  The  computational  procedure  applied  is  as  out- 
lined in  previous  S^  reports  (e.g.,  Bache,  et  al.,  1975,  1976; 
Barker,  et  al. , 1976).  The  methods  have  been  shown  to  give 
theoretical  seismograms  that  agree  quite  well  with  observa- 
tions when  the  appropriate  model  parameters  are  included  in 
the  calculations. 

The  coupling  of  the  explosion-produced  shock  waves  into 
elastic  waves  is  represented  by  a reduced  displacement  poten- 
tial which  provides  an  equivalent  elastic  source.  A brief 
description  of  the  constitutive  model  and  material  properties 
for  the  explosions  in  salt  is  given  in  Section  II.  Four 
calculations  were  carried  out  in  which  all  material  proper- 
ties were  identical  except  for  the  overburden  pressure  which 
varies  proportional  to  the  burial  depth.  We  find  that  the 
long  period  level  of  the  source  spectrum  is  rather  insensitive 
to  overburden  pressure.  However,  at  the  low  confining  pressures 
characteristic  of  shallow  burial  depths,  tension  cracks  open 
in  the  material  and  the  source  spectrum  is  strongly  peaked 
near  1 Hz,  the  range  of  interest  for  body  wave  magnitude. 

In  Section  III  we  present  our  theoretical  analysis  of 
the  nij^-yield  relationships  for  explosions  in  salt.  We  begin 
by  computing  theoretical  seismograms  for  a series  of  explo- 
sions varying  in  yield  and  burial  depth.  The  source-receiver 
travel  path  is  meant  to  be  appropriate  for  explosions  in 
Eurasia.  We  find  an  m^-log  yield  relationship  with  slope 
of  2U30ut  unity,  though  it  can  vary  from  this  value  depending 
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on  the  population  being  fit.  Further,  the  theoretical 
values  agree  rather  well  with  observations. 

Our  next  step  was  to  compute  theoretical  seismo- 
grams for  several  hypothetical  salt  events  at  NTS.  For 
NTS  events  we  have  estimates  for  the  source-receiver 
path  characteristics  and  can  accurately  predict  m^^.  The 
major  difference  between  the  NTS  emd  Eurasian  path  is  in 
the  attenuation  which  is  characterized  by  a parameter  called 
t*.  In  this  case  the  t*  values  are  1.06  and  0.6  for  the  two 
paths.  The  nij^-log  yield  curve  for  the  hypothetical  salt 
events  is  then  0.5-0. 6 m^^  units  below  that  for  otherwise 
identical  events  in  Eurasia.  We  caution  that  this  dif- 
ference is  based  on  our  "first-cut"  estimate  for  the  source- 
receiver  path  for  Eurasian  events.  But  while  we  are  not 
sure  exactly  what  the  irij^  difference  between  identical  events 
at  the  two  sites  is,  we  believe  it  is  substantial. 


The  next  step  in  our  analysis  is  to  compare  the  nij^- 
yield  data  for  the  salt  events  to  that  for  events  in  granite 
at  NTS.  We  find  that  the  granite  itij^-log  yield  curve  lies 
midway  between  those  for  the  Eurasian  and  NTS  salt  events. 
This  reflects  the  fact  that  the  granite  is  a somewhat  higher 
coupling  material  than  the  salt. 

Finally,  we  look  at  the  event  SALMON,  a U.S.  explosion 
in  salt.  The  comparison  of  theoretical  and  actual  m^  is 
somewhat  eunbiguous  because  of  the  scatter  in  the  observed 
values.  However,  we  find  that  with  t*  = 0,6  our  theoretical 
seismogram  at  the  best  of  the  teleseismic  stations  (COL) 
agrees  quite  well  with  the  observations.  On  the  other 
hand,  higher  t*  values  give  better  agreement  with  the 
average  m.  over  the  network  based  on  many  observations. 
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II.  THE  EQUIVALENT  ELASTIC  SOURCE  FOR  EXPLOSIONS  IN  SALT 

Coupling  of  the  explosion-produced  shock  wave  was  com- 
puted with  the  one-dimensional  (spherically  symmetric)  finite 
difference  code,  SKIPPER.  Description  of  the  technique  and 
the  constitutive  models  may  be  found  in  Cherry,  et  al.  (1975) , 

The  material  parameters  for  the  calculation  were: 

P wave  velocity,  a = 4.2  km/sec, 

S wave  velocity,  B = 2.41  km/ sec. 

Density,  p =2.15  gm/cm*. 

Porosity,  0 =0.0, 

Failure  envelope  parameters: 

Y = 0.066  kbar, 

0 

= 0.467  kbar, 

= 0.42  kbar. 


m 


m 


The  failure  envelope  parameters  satisfy 


Y(P)  = Y + Y 


m 


m ' m ' 


' f ' 


+ X. 


m' 


P > P, 


m 


(2.1) 


P - 


where  J"  is  the  third  deviatoric  stress  invariant  and  P 

3 

is  the  pressure  including  the  overburden.  The  Y is  the 
maximum  stress  difference  or  twice  the  maximum  shear  stress. 
The  meeming  of  these  parameters  is  indicated  schematically 
in  Figure  2.1. 

These  strength  parameters  give  an  unconfined  compres- 
sive strength  of  0.411  kb  and  a maximum  allowable  stress  dif- 
ference of  0.533  kb. 
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The  finite  difference  calculation  of  the  outgoing 
shock  wave  is  carried  into  the  elastic  regime.  From  the 
elastic  field  an  equivalent  elastic  point  source  is  com- 
puted in  the  form  of  the  reduced  displacement  potential  (4') 
which  is  related  to  the  displacement  by 


u(R,t) 


y(T)  . jw 

r2  Ra 


(2.2) 


T « t - R/a  . 

A convenient  display  of  the  elastic  source  is  given  by  the 
amplitude  of  the  Fourier  transformed  reduced  velocity  poten- 
tial, l’?(w)|. 

Four  calculations  were  carried  out  with  only  the  over- 
burden pressure  varying.  Overburden  pressure  is  simply  pgH 
where  H is  the  explosion  burial  depth  and  g is  the  ac- 
celeration of  gravity.  The  source  functions  for  these  four 
calculations  are  shown  in  Figure  2.2.  Note  that  the  depth 
varies  from  quite  shallow  (0.15  km)  to  quite  deep  (2.5  km). 

In  comparing  the  source  functions  in  Figure  2.1 

there  are  two  frequency  ranges  of  primary  interest.  The 

teleseismic  surface  wave  magnitude  (M  ) depends  on  the 

s 

amplitude  of  the  source  function  around  0.05  Hz.  On  the 
other  hand,  for  irij^  the  imgortant  frequencies  are  around 
1 Hz.  The  value  of  the  |'i'(a))|  at  these  two  frequencies  is 
tabulated  in  Table  2.1. 

From  the  tabulated  values  we  see  that  the  long 
period  level  is  mildly  sensitive  to  depth,  varying  by  only 
20  percent  over  the  entire  range.  However,  the  short  period 
level  varies  by  a factor  of  2.2  from  largest  to  smallest. 

This  is  because  of  the  marked  peaking  of  the  source  function 
at  shallow  depths.  This  peaking  is  due  to  the  opening  of 
tensile  cracks  in  the  material  at  low  confining  pressures. 


0.15  (251)/ 


= 0 . 5_J^5,5-)-^  \ 
=-^3  ■(issT'^  ^\\ 

2.5  (252K  A 
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Figure  2.2. 


Frequency  (Hz) 


Source  functions  for  salt  with  variable  over- 
burden (indicated  by  the  depth  H,  in  km) . Th 
frequency  axis  is  scaled  to  150  kt,  the  ampli 
tude  axis  to  0.02  kt. 
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In  shorty  the  body  wave  portion  of  the  spectrum  is  much 
more  sensitive  to  depth  of  burial  than  the  long  period  or 
surface  wave  portion. 

TABLE  2.1 

SOURCE  FUNCTION  AMPLITUDE  AT  CHARACTERISTIC  BODY 


AND  SURFACE  WAVE 

FREQUENCIES 

Source 

Depth  (km) 

A 

V (.05  Hz,mM 

(1  Hz 

251 

0.15 

8.2 

15.1 

255 

0.5 

7.9 

8.5 

253 

1.3 

7.4 

7.4 

252 

2.5 

6.8 

6.8 

The  value  of  4'  («)  obtained  from  free  field  measure- 
ments for  SALMON  (5  kt,  823  m DOB)  was  approximately  3100  m 
(Patterson,  1966).  This  ¥ {<»)  cube  root  scales  to  12.4  m^ 
for  a yield  of  0.02  kt.  From  Table  2.1  the  comparable 
computed  values  (4'  (0.05  Hz))  is  a factor  of  1.6  below  that 
for  SALMON.  This  corresponds  to  a magnitude  difference 
between  the  SALMON  and  calculated  sources  of  approximately 
0.2  magnitude  units. 

No  material  strength  data  were  available  for  either 
the  GNOME  or  SALMON  salt.  The  strength  parameters  used  for 
the  source  calculations  were  obtained  from  Heard,  et  al. 
(1975)  and  Lyuke,  et  al.  (1976).  The  Lawrence  Livermore 
Laboratory  data  were  obtained  from  an  artificial  aggregate 
of  hydrostatically  compressed  polycrystalline  salt  with  a 
purity  of  99.4  to  99.9  percent.  The  calculation  of  the 
SALMON  equivalent  source  should  be  performed  with  strength 
data  appropriate  to  the  SALMON  salt. 


I 

"i 

i III.  THE  DEPENDENCE  OF  Itij^  ON  YIELD  FOR  EXPLOSIONS  IN  SALT 

I 3.1  COMPUTATIONAL  METHOD 

I 

I The  computational  method  for  computing  body  wave 

^ j seismograms  includes  the  following  elements: 

1.  The  reduced  displacement  potential  (RDP)  is 

* computed  by  our  one-dimensional  finite  dif- 

ference code  (SKIPPER) . The  results  were  pre- 
sented in  Section  II. 

t 

2.  The  RDP  is  input  to  a prograun  which  computes 
the  detailed  crustal  reverberations  for  an 

; i equivalent  elastic  point  source  at  depth  in 

i I a plane-layered  crustal  model  (Bache  and 

Harkrider,  1976] . 

3.  Crustal  reverberations  in  a plane-layered  model 
of  the  crust  at  the  receiver  are  computed  using 
the  method  of  Haskell  [1962] . 

4.  The  effect  of  the  upper  mantle  is  computed  by 
generalized  ray  theory  as  implemented  by 
Helmberger  and  described  in  Wiggins  and 
Helmberger  [1974] . 

5.  Anelastic  attenuation  and  the  attendant  disper- 
sion are  included  via  the  operator 

exp  - irf  t*  |l  - I i An  (^3^)  , 

using  the  formulation  suggested  by  Strick  [1970] . 
Here  f is  frequency  and  the  controlling  para- 


6. 


meter  is  t*  = T/Q,  the  ratio  of  travel  time  to 
the  average  path  material  quality  factor. 


The  response  of  the  appropriate  instrument  is 
convolved  with  the  ground  motion. 
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3.2  EARTH  MODELS 

Two  models  for  the  source  crustal  structure  were  used. 
These  are  tabulated  in  Table  3.1.  For  the  crust  in  the 
receiver  region  a simple  model  having  no  important  effect  on 
the  frequency  content  of  the  seismogram  was  chosen.  However, 
the  absolute  amplitude  is  roughly  proportional  to  the  inverse 
of  the  velocity  of  the  top  layer  (see  Appendix  B of  Barker, 
et  al.  [1976]).  This  model  is  tabulated  in  Table  3.2.  The 
crustal  models  at  the  source  and  receiver  accounted  for  the 
top  20  km  of  the  earth.  We  also  need  a model  for  the  crust 
and  upper  mantle  below  this  depth. 

The  stations  of  interest  for  the  synthetic  seismo- 
gram calculations  were  at  great  distances  from  the  source 
(>  45®).  At  this  distance  the  upper  mantle  does  little  more 
than  geometrically  spread  the  single  ray  that  makes  up  the 
P wave  arrival.  The  amplitude  is  controlled  by  the  struc- 
ture near  the  tiarning  point  which  is  at  great  depth  (>  1000 
km) . In  selecting  an  upper  mcintle  model  then  there  are  two 
points  to  be  kept  in  mind.  First,  there  is  little  evidence 
of  lateral  inhomogeneity  in  the  mantle  at  these  depths. 

Second,  most  proposed  earth  models  predict  about  the  seime 
amplitude  (geometric  spreading)  for  the  distances  of  interest. 

In  view  of  the  comments  of  the  previous  paragraph,  the 
upper  mantle  model  selected  for  our  calculations  was  one  of 
the  best  available  average  earth  models.  This  model  is  C2 
proposed  by  Anderson  and  Hart  [1976] , which  we  refer  to  as 
C2AH. 

The  remaining  earth  model  quantity  to  be  selected  is 
t*.  For  these  computations  we  use  t*  = 0.6.  We  believe  that 
t*  for  Eurasian  travel  paths  is  considerably  lower  than  for 
paths  such  as  NTS-Alaska  that  we  have  extensively  studied 
(e.g.,  Bache,  et  al. , 1977).  For  the  latter  paths  t*  approxi- 
mately equal  to  1.0  is  appropriate.  The  value  of  0.6  is  an 
estimate  that  requires  more  research  for  refinement. 
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TABLE  3.1 

SOURCE  REGION  CRUSTAL  MODELS* 


Depth  (km) 

a (km/sec) 

B (km/sec) 

p (qm/cm 

Model  SSI 

0.1 

1.4 

0.8 

1.6 

3.3 

4.2 

2.4 

2.15 

o 

. 

in 

4.6 

2.6 

2.5 

00 

. 

o 

5.0 

2.7 

2.7 

— 

6.0 

3.5 

2.8 

Model  SS2 

0.18 

1.4 

0.8  . 

1.6 

0.24 

' 5.0 

2.7 

2.6 

3.3 

4.2 

2.4 

2.15 

5.0 

4.6 

2.6 

2.5 

00 

. 

o 

5.0 

2.7 

2.7 

— 

6.0 

3.5 

2.8 

I 

i 


» 

< 

I 


} 


f 


These  structvjres  are  not  meant  to  represent  any 
particular  area  but  merely  to  be  typical  of  what 
is  encountered  in  salt  dome  regions  (e.g.,  Warren, 
et  al. , 1966) . 


TABLE  3.2 

RECEIVER  REGION  CRUSTAL  MODEL 


Depth  (km) 

a (km/sec) 

8 (km/sec) 

p (qm/cm*) 

2.58 

3.67 

2.31 

2.40 

4.84 

5.42 

3.27 

2.60 

11.61 

5.80 

3.45 

2.60 

— 

6.00 

3.50 

2.80  \ 

S 
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3.3  SYNTHETIC  SEISMOGRAMS 

Synthetic  seismograms  were  computed  at  two  epicentral 
distances,  56.6*  and  66.8*.  The  peurameters  for  the  calcula- 
tion axe  summarized  in  Table  3.3. 

Note  that  the  source  crustal  structure  denoted  SS2 
in  Tcible  3.1  was  used  for  all  the  synthetic  seismograms 
except  the  two  at  the  smallest  burial  depth.  For  these  two 
the  high  velocity  cap  on  the  salt  layer  was  removed  to  make 
up  the  structure  SSI. 

We  must  say  a few  words  about  the  P-pP  lag  time  in 
the  calculation  and  its  relation  to  the  burial  depth.  The 
P and  S wave  velocities  specifying  the  structures  in  Tables 
3.1  and  3.2  are  elastic  wave  velocities  consistent  with  well 
log  or  laboratory  measurements.  It  is  well-known  that  the 
initial  ground  motion  at  ground  zero  almost  always  arrives 
later  than  would  be  expected  from  the  sonic  velocity  logs 
(the  energy  travels  part  of  the  way  as  a shock  wave  with 
velocity  higher  than  the  sound  speed) . Also,  study  of  the 
P-pP  lag  from  far-field  P wave  recordings  indicates  a greater 
lag  than  expected  from  the  elastic  analysis  (e.g..  Springer, 
1974) . In  the  elastic  analysis  the  P-pP  lag  may  be  computed 
from 


?-pP  lag  - 2^  ^ cos  , 
i-1  ^ 


9 . = sin" 


(3.1) 


where  is  the  thickness  of  the  layer  of  velocity  be- 
tween the  source  and  the  free  surface.  The  ray  parameter, 
p,  is  determined  by  the  generalized  ray  theory  calculations. 
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TABLE  3.3 

PARAMETERS  FOR  SYNTHETIC  SEISMOGRAM  CALCULATIONS 


Yield 

Depth* 

P-pP  lag* 

Source 

Crustal 

Structure 

1 

0.165 

0.19 

251  (0.15) 

SSI 

8 

0.165 

0.19 

251  (0.5) 

SSI 

25 

0.60 

0.50 

255  (0.5) 

SS2 

65 

0.60 

0.50 

255  (0.5) 

SS2 

15 

1.15 

0.80 

253  (1.3) 

SS2 

35 

1.15 

0.80 

253  (1.3) 

SS2 

40 

2.4 

1.49 

252  (2.5) 

SS2 

65 

2.4 

1.49 

252  (2.5) 

SS2 

200 

1.0 

0.72 

253  (1.3) 

SS2 

200 

1.5 

1.00 

253  (1.3) 

SS2 

200 

2.0 

1.27 

252  (2.5) 

SS2 

200 

2.5 

1.55 

252  (2.5) 

SS2 

The  depth  used  in  the  calculations  was  modified  to  increase 
the  P-pP  lag,  see  the  text. 

** 

X The  source  with  overburden  depth  closest  to  the  shot  depth 

was  selected  since  calculations  were  not  done  at  all  source 
depths . 
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In  this  case  p = 0.06  sec/km.  This  corresponds  to  0 = 14.7® 
in  the  source  layer  (a  = 4.2  km/sec) . 

It  is  clear  from  (3.1)  that  to  increase  the  P-pP  lag 
we  must  either  increase  or  decrease  a^.  The  latter 
chemge  will  affect  the  reflection  coefficients  as  well  as 
the  lag.  Therefore,  we  chose  to  increase  the  to  obtain 
lags  that  are  somewhat  greater  than  those  predicted  by  the 
given  elastic  velocities  and  depths.  For  this  reason  the 
actual  calculations  were  done  with  the  depth  in  Table  3.3 
multiplied  by  1.2. 

The  computed  synthetic  seismograms  are  shown  in 
Figures  3.1  and  3.2.  The  source  ftinction,  yield  and  depth 
for  each  calculation  is  printed  on  the  records.  Also  ap- 
pearing on  the  records  are  the  values  from  measurements 
of  the  b and  A phases.  The  b phase  is  the  first  cycle  on 
the  record.  The  A phase  is  the  meucimum  phase  in  the  first 
few  cycles  and  is  indicated  with  a bar.  The  is  computed 
from 

mj^  - log  I + B,  (3.2) 

where  B is  taken  to  be  3.5  for  the  station  at  56.6®  and 
3.7  for  the  station  at  66.8®.  The  eunplitude  in  (3.2)  is 
measured  peak-to-peak.  Given  in  parentheses  with  each  m^^ 
is  the  period  of  the  cycle  measured.  The  period  is  seen 
to  be  dependent  on  the  explosion  depth  and  yield.  The  raw 
data  for  the  seismogreuns  of  Figures  3.1  and  3.2  is  summarized 
in  Tcdsle  3.4. 

3.4  THE  DEPENDENCE  OF  m^  ON  YIELD 

The  next  step  is  to  examine  the  dependence  of  on 
yield.  For  this  purpose  we  plot  mj^  versus  log  W for  m^ 
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Figure  3.1.  Synthetic  seismogrzuns  at  56.6*.  The  numbers  at 
the  left  are  ground  motion  in  microns  at  1 Hz. 
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TABLE  3.4 


AMPLITUDE/PERIOD 

DATA  FOR  THE 

SYNTHETIC 

SEISMOGRAMS 

OF 

FIGURES  3.1 

AND  3.2 

Yield 

P-pP  Lag 

b 

!b 

A 

!a 

Range  » 

56. 6«  (B  = 

3.5) 

1 

0.19 

0.85 

0.52 

1.27 

0.61 

8 

0.19 

8.32 

0.56 

12.0 

0.58 

25 

0.50 

38.7 

0.64 

103.0 

0.83 

65 

0.50 

107.0 

0.65 

266.0 

0.82 

15 

0.80 

20.1 

0.68 

34.1 

0.81 

35 

0.80 

46.6 

0.68 

81.5 

0.82 

40 

1.49 

44.7 

0.65 

69.4 

0.78 

65 

1.49 

74.1 

0.66 

117.0 

0.79 

200 

0.72 

308.0 

0.72 

632.0 

0.89 

200 

1.00 

311.0 

0.76 

293.0 

0.73 

200 

1.27 

237.0 

0.68 

351.0 

0.78 

200 

1.55 

230.0 

0.67 

362.0 

0.79 

Range  « 

66.8-  (B  - 

3.7) 

1 

0.19 

0.75 

0.52 

1.07 

0.59 

8 

0.19 

7.36 

0.56 

10.6 

0.56 

25 

0.50 

33.3 

0.63 

86.8 

0.83 

65 

0.50 

91.9 

0.64 

224.0 

0.82 

15 

0.80 

16.9 

0.67 

26.1 

0.76 

35 

0.80 

40.6 

0.68 

55.7 

0.73 

40 

1.49 

38.4 

0.64 

52.6 

0.75 

65 

1.49 

63.7 

0.65 

88.7 

0.76 

200 

0.72 

267.0 

0.72 

474.0 

0.86 

200 

1.00 

269.0 

0.76 

269.0 

0.76 

200 

1.27 

206.0 

0.68 

284.0 

0.77 

200 

1.55 

196.0 

0.66 

295.0 

0.78 
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2md  njj^  at  the  two  ranges.  The  best  least  squares  linear 
fit  to  the  data  was  computed  giving  relations  of  the  form 

= a^  log  W + a^,  (3.3) 

The  regression  was  done  assuming  no  error  in  the  yield 
values.  The  results  are  summarized  in  Table  3.5.  In 
Figures  3. 3-3. 4 the  m^  is  plotted  for  the  two  stations. 

The  nij^-log  yield  values  are  separated  into  two  groups 
according  to  their  scaled  depth  of  burial.  The  shallow 
events  have  scaled  depths  DG  [82.5,  205],  where  the  actual 
depth  H a The  deeply  buried  events  have  D varying 

between  351  and  701  meters. 

From  the  graphs  and  Table  3.5,  we  see  that  there  is 
almost  no  difference  between  the  seismogreuns  at  the  two 
ranges.  At  these  ranges  the  upper  mantle  spreading  is  nearly 
independent  of  frequency  and  proportional  to  A”*.  This 
spreading  gives  log  A/T  values  at  the  farther  station  that 
are  about  0.07  smaller  than  those  at  the  nearer  station. 
However,  the  distance  correction  (B)  differs  by  0.2  magni- 
tude units.  That  is,  the  amplitude  decay  predicted  by  the 
upper  mantle  model  is  not  consistent  with  the  variation  of 
B in  this  distance  range.  Thus,  the  magnitude  at  the  farther 
station  should  be  about  0.13  magnitude  units  larger  than  that 
at  the  nearer  and  this  is  consistent  with  the  values  we  get. 

3.5  COMPARISON  TO  OBSERVATIONS 

Several  of  the  theoretical  events  are  similar  in  depth 
and  yield  to  actual  events.  The  m^^  for  these  events  can  be 
compared  to  the  theoretical  values  plotted  in  Figures  3.3 
and  3.4.  The  agreement  between  the  theoretical  and  actual 
is  generally  quite  good. 

3.6  COMPARISON  TO  HYPOTHETICAL  SALT  EVENTS  AT  NTS 

We  have  a we 11 -developed  capability  to  compute  theo- 
retical seismograms  for  NTS  events  that  have  m^  quite  close 
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TABLE  3 . 5 

! 

BEST  LINEAR  RELATIONS  FOR  THE  DATA  OF  FIGURES 
i 3.1  AND  3.2 

Data  Set 


i A = 56.6* 

f 

1 Entireset 

[ 

= 1.02  log  W + 3.95 

* 1.01  log  W + 3.93 

I Da82.5,205] 

= 1.07  log  W + 3.93 

= 1.14  log  W + 3.86 

Da351,7011 

it^  = 0.95  log  W + 4.05 

= 0.92  log  W + 4.03 

I A » 66.8* 

Entire  set 

= 1.02  log  W + 3.90 

» 1.00  log  W + 4.06 

Da 82. 5, 205] 

= 1.07  log  W + 3.88 

it^  = 1.12  log  W + 4.01 

Da351,701] 

= 0.95  log  W + 3.99 

» 0.92  log  W + 4.12 

1 
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i 
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1 
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1 
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yield  at  56.6*.  Least  square  lines  are  fit  to  the  total  data 
I two  subgroups  divided  by  scaled  depth  of  burial. 
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to  that  observed.  Therefore,  let  us  construct  theoretical 
seismograms  for  events  just  like  those  discussed  in  the  pre- 
vious section,  changing  the  path  parameters  to  be  appropriate 
for  NTS  events.  Once  again  the  upper  mantle  model  C2AH  is 
used,  but  the  range  is  33.4*.  The  t*  = 1.06.  With  these 
path  parameters  a B value  of  3.32  gives  theoretical  m^^  that 
closely  agree  with  network  averages  for  NTS  events. 


Theoretical  seismograms  for  hypothetical  NTS  salt 
events  were  computed  for  five  of  the  twelve  events  described 
in  Table  3.3.  The  seismograms  are  shown  in  Figure  3.5.  The 
data  is  simmarized  in  Table  3.6.  These  data  should  be  com- 
pared to  the  corresponding  data  in  Table  3.4. 


The  nij^  for  the  NTS  travel  path  events  are  compared 
to  those  for  the  Eurasian  travel  path  events  in  Table  3.7. 

The  comparison  is  given  for  both  the  itij^  and  m^  values.  The 
difference  is  computed  using  the  average  for  the  two  Eurasian 
travel  paths.  Except  for  the  shallow  low  yield  event,  the 
nij^  difference  is  fairly  constant,  varying  by  less  than  0.1. 
For  m^  the  15  kt  event  is  also  anomalous  because  the  apparent 
period  of  the  A phase  is  quite  large  due  to  the  smooth  inter- 
ference between  P and  pP  at  the  low  t*.  The  instrument 
response  correction  is  based  on  this  large  apparent  period 

A 

and  causes  the  for  the  NTS  path  to  be  artifically  large. 


3.7  COMPARISON  TO  EVENTS  IN  GRANITE 

It  is  interesting  to  compare  nij^-yield  for  the  salt 
events  to  that  for  explosions  in  granite.  The  representative 
event  in  granite  is  PILEDRIVER.  We  have  computed  synthetic 

seismograms  for  PILEDRIVER  in  the  past  (e.g.,  Bache,  et  al. 

A ' 

[1975])  and  obtain  a theoretical  m^^  of  5.55.  The  calcula- 
tion used  t*  = 1.06  and  the  S2une  NTS  - 33.4®  path  as  the 
seismograms  in  the  previous  section. 

Let  us  first  point  out  the  major  factors  influencing 
the  theoretical  nij^  for  PILEDRIVER.  Our  one-dimensional 
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Figure  3.5.  Theoretical  seismogr2uns  for  hypothetical  events 
in  salt  at  NTS.  The  source  parameters  are  the 
same  as  for  the  corresponding  Eurasian  salt 
events  studied  in  previous  sections.  The  path 
parameters  are  appropriate  for  an  NTS-Alaska 
path. 


22 


R-3057 


w 


C 


TABLE  3.6 

AMPLITUDE/PERIOD  DATA  FOR  THE  SYNTHETIC  SEISMOGRAMS 

OF  FIGURE  3.5 


Yield 

Range  » 

P-pP  Lag 

33.4*  (B  « 

b 

3.32) 

!b 

A 

1 

0.19 

0.32 

0.69 

0.59 

0.79 

25 

0.50 

24.6 

0.81 

48.1 

0.87 

15 

0.80 

16.2 

0.89 

47.8 

1.17 

40 

1.49 

30.8 

0.84 

42.2 

0.92 

200 

1.27 

157.0 

0.84 

173.0 

0.85 
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TABLE  3.7 

COMPARISON  OF  FOR  NTS  AND  EURASIAN  TRAVEL  PATHS 

0 


Eurasian 

Paths 

NTS  Path 

Yield 

^=56.6" 

A=66.8® 

A=33.4“ 

Difference 

A 

“b 

1 

3.82 

3.96 

3.20 

0.69 

25 

5.59 

5.72 

5.06 

0.59 

15 

5.12 

5.24 

4.98 

0.20 

40 

5.45 

5.55 

4.98 

0.52 

200 

6.15 

6.27 

5.63 

0.58 

b 

1 

3.91 

3.86 

2.99 

0.90 

25 

5.48 

5.42 

4.81 

0.64 

15 

5.17 

5.10 

4.58 

0.56 

40 

5.54 

5.48 

4.89 

0.62 

200 

6.24 

6.18 

5.59 

0.62 

24 
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source  calculation  for  PILEDRIVER  gives  a source  function 

that  is  larger  than  that  for  salt.  This  is  illustrated  in 

Figure  3.6  where  the  largest  (251)  and  smallest  (252)  of  the 

salt  source  functions  are  compeured  to  that  for  PILEDRIVER 

granite.  All  other  things  being  equal,  we  expect  the  ampli- 

tude  of  the  telsesismic  body  waves  to  scale  with  a ¥ where 

«e  ® 

a xs  the  source  region  P wave  velocity  and  ’F  is  the  ampli- 
s 

tude  of  the  source  spectrvan  at  the  appropriate  frequency, 

about  1 Hz.  For  salt,  a =4.2  km/sec  and  for  granite  a = 

s s 

5.3  km/sec.  From  this  first-order  scaling  we  would  expect 
the  granite  body  wave  amplitude  to  be  a factor  of  2-3  higher 
them  that  for  otherwise  identical  events  in  salt. 

The  pronounced  peak  in  the  PILEDRIVER  source  func- 
tion is  clearly  responsible  for  much  of  the  coupling  dif- 
ference between  the  PILEDRIVER  granite  and  the  salt.  This 
peaking  is  caused  by  several  peculiarities  of  the  PILEDRIVER 
granite  — particularly  the  presence  of  water  and  pre-exist- 
ing cracks  in  the  granite  (see  Bache,  ^ al. , 1975).  These 
conditions  need  not  apply  at  other  granite  sites  or  at 
greater  depths  at  the  same  site. 

A second  factor  strongly  controls  the  PILEDRIVER  nij^. 
This  event  is  characterized  by  an  exceptionally  small  lag 
between  P and  pP.  This  causes  destructive  interference 
and  a smaller  m^^  than  expected  from  source  coupling  alone. 

In  order  to  construct  an  itij^-yield  curve  for  granite 
similar  to  those  shown  for  salt,  we  need  several  more 
granite  source  functions  that  take  account  of  the  depth  ef- 
fects on  source  coupling.  While  these  calcualtions  are  in 
progress,  we  can  anticipate  pretty  closely  what  the  results 
will  be.  As  we  increase  the  depth,  the  source  function  would 
get  smaller  because  of  the  increasing  overburden  pressure. 

We  also  might  reasonably  reduce  the  amount  of  precracking  in 
the  material.  In  any  case,  the  major  effect  would  be  on 
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Figure  3.6«  The  PILEDRIVER  source  function  compared  to  the 
largest  and  smallest  of  the  salt  source  func- 
tions from  Figure  2.2.  The  frequency  axis  is 
scaled  to  150  kt,  the  amplitude  eocis  to  0.02  kt 
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the  1 Hz  portion  of  the  spectrum  which  would  be  substantially 
reduced  with  increasing  depth  as  was  the  case  for  salt 
(Figure  2.1).  At  the  same  time  the  amount  of  destructive 
interference  between  P and  pP  becomes  less  with  depth.  The 
result  should  be  an  nij^  log  W curve  that  has  slope  of  nearly 
unity. 

The  nij^-yield  values  for  salt-Eurasiam  travel  path, 
salt-NTS  travel  path  and  PILEDRIVER  granite-NTS  travel  path 
are  compared  in  Figure  3.7.  Three  lines  of  unit  slope  are 
shown  through  the  data.  The  upper  most  line  is  the  least 
squares  fit  to  the  twelve  salt-Eurasian  path  events  plotted 
in  Figures  3.3  and  3.4.  The  m^^  values  are  the  average 
values  from  the  two  stations.  The  lowest  line  is  0.55 
magnitude  vinits  below  the  salt-Eurasian  path  line.  This  is 
representative  of  the  difference  between  the  salt-Eurasian 
path  and  salt-NTS  path  nij^  values  compared  in  Table  3.7. 

The  unit  slope  line  through  the  PILEDRIVER  theoreti- 
cal falls  almost  equidistant  between  the  lines  for  the 

two  populations  of  salt  events.  That  is  to  say,  events  like 
PILEDRIVER  are  expected  to  give  m^^  values  that  are,  on  the 
average,  0.1-0. 3 magnitude  units  above  those  for  similar 
NTS  events  in  salt.  On  the  other  hand,  the  PILEDRXVER-like 
events  give  values  that  are  about  0.1-0. 3 units  lower 
than  Eurasian  salt  events. 

3.8  SALMON 

There  has  been  one  U.S.  explosion  in  salt  that  was 
fairly  well  recorded  teleseismically.  This  is  the  ::  5 kt 
explosion  SALMON  in  Mississippi.  While  discussing  salt  ex- 
plosions we  should  compare  our  theoretical  seismograms  to 
those  for  SALMON. 

Synthetic  seismograms  were  computed  for  SALMON  at 
A > 48.8*.  The  source  region  crustal  structxire  is  tabulated 
in  Table  3.8.  The  receiver  region  structure  is  the  same  as 
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that  used  in  all  calculations  in  this  report,  that  specified 
in  Tcible  3.2. 

The  actual  depth  of  SALMON  was  0.83  km.  This  depth 
is  midway  between  those  for  the  salt  source  functions  253 
and  255.  We  chose  to  compute  seismograms  using  255,  the 
larger  of  the  two.  For  the  actual  seismogram  calculations 
we  again  increased  the  depth,  this  time  by  28  percent  to 
1.06  km.  The  result  is  a P-pP  delay  of  0.56  sec.  This  is 
neeurly  the  seune  as  the  delay  time  of  0.58  seconds  estimated 
by  Springer  [1974]  from  the  arrival  time  at  ground  zero. 

The  earth  model  for  the  upper  mantle  path  was  again 
C2AH,  this  time  with  A » 48.8®  (P  = 0.072  secAm)  , the  range 
from  SALMON  to  the  WWSSN  station  at  College,  Alaska  (COL) . 

As  was  the  case  with  the  Eurasian  paths,  it  is  difficult  to 
confidently  select  t*  for  this  unfamiliar  travel  path.  We 
suspect  that  the  appropriate  value  is  smaller  than  for  NTS 
events.  Seismograms  were  computed  for  several  values  of  t*. 
The  standard  WWSSN  seismograph  response  was  convolved  with 
the  ground  motion. 

The  theoretical  seismograms  for  SALMON  are  shown  in 
Figure  3.8.  The  data  from  the  seismograms  is  summarized  in 
Table  3.9.  For  the  calculations  B = 3.5. 

How  do  the  theoretical  seismograms  compare  with  obser- 
vations? The  teleseismic  observations  of  SALMON  cure  sxim- 
marized  by  Jordan,  et  al.  [1966] . The  magnitude  values 
given  by  these  authors  are  shown  plotted  versus  distemce 
in  Figure  3.9.  Restricting  attention  to  those  values  beyond 
16®,  the  average  magnitude  is  found  to  be  m^^  = 4.35.  How- 
ever, most  of  the  teleseismic  magnitudes  (A  » 30®)  fall  above 
this  average  value. 

Jordan,  et  al.  show  seismograms  for  a number  of  tele- 
seismic stations.  The  best  recording  is  that  at  COL,  the 
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Figure  3.8.  Synthetic  seismograms  for  SALMON  at  A = 48.8 
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station  for  which  our  synthetics  were  computed.  At  this 
station  A = 20.5,  » 0.8  and  m^  » 4.97.*  Comparing  to 

the  data  of  Table  3.9,  we  see  that  the  theoretical  seismo- 
gram with  T/Q  = 0.6  is  in  excellent  agreement  with  the  ob- 
servation. On  the  other  hand,  higher  values  of  t*  give 
better  agreement  with  the  average  m^^  for  the  event. 


Jordam,  ^ al.  use  B » 3.7  at  this  distance,  a value  not 
in  agreement  with  the  LRSM  shot  report  tables. 


( 
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